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SGd,b{n) with d = 2 and & = 3, 4 are also obtained. 

Keywords: Spanning trees, Sierpinski gasket, exact solutions. 



I. INTRODUCTION 

The enumeration of the number of spanning trees Nst^G) on a graph 
G was first considered by Kirchhoff in the analysis of electric circuits 
It is a problem of fundamental interest in mathematics H, Q |^ and 
physics 0, Q ■ The number of spanning trees is closely related to thepar- 
tition function of the g-state Potts model in statistical mechanics l^- 
Some recent studies on the enumeration of spanning trees and the calcula- 
tion of their asymptotic grow th constants on regular lattices were carried 
out in Refs. 0, 0, IT^ Il3l|. It is of interest to consider spanning trees 
on self-similar fractal lattices which have scaling invariance rather than 
translational invariance. Fractals are geometric structures of noninteger 
Hausdorff dimension realized by repeated construction of an elementary 
shape on progressively smaller length scales jQ, . A well-known exam- 
ple of fractal is the Sierpinski gasket which has been extensively studied 
in several contexts d iHJl, 19 , 20, 21, 22, 23., 2± „25,, ,26]. We shaU 
derive rigorously the numbers of spanning trees on the Sierpinski gasket 
with dimension equal to two, three and four. The corresponding asymp- 
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totic growth constants have simple expressions. We shall also conjecture 
the general expression of the number of spanning trees on the Sierpinski 
gasket with arbitrary dimension. 



II. PRELIMINARIES 



We first recall some relevant definitions for spanning trees and the Sier- 
pinski gasket in this section. A connected graph (without loops) G = (V, E) 
is defined by its vertex (site) and edge (bond) sets V and E '2I, [23| ■ Let 
v{G) = \V\ be the number of vertices and e{G) = \E\ the number of edges 
in G. A spanning subgraph G" is a subgraph of G with the same vertex 
set V and an edge set E' C E. As a tree is a connected graph with no 
circuits, a spanning tree on G is a spanning subgraph of G that is a tree 
and hence e(G") = v(G) — 1. The degree or coordination number ki of a 
vertex Vi G V is the number of edges attached to it. A /c-regular graph is 
a graph with the property that each of its vertices has the same degree k. 
In general, one can associate an edge weight Xa to each edge connecting 
adjacent vertices Vi and Vj (see, for example |liJ|). For simplicity, all edge 
weights are set to one throughout this paper. A well-known method to enu- 
merate spanning trees is to construct the Laplacian matrix Q{G), which is 
the v{G) X v{G) matrix with element Q{G)ij — kiSij — A{G)ij where A{G) 
is the adjacency matrix. One of the eigenvalues of Q{G) is always zero. 
Denote the rest as A(G)i, 1 < i < v{G) — 1, then the number of spanning 

trees Nst{G) = {1/viG)) UtT' HG)r i- 

When the number of spanning trees Nst{G) grows exponentially with 
v(G) as v{G) — *■ cxD, there exists a constant zq describing this exponential 
growth HQ: 

ZG = lim — (2.1 

i>(G)^oo V(G) 

where G, when used as a subscript in this manner, implicitly refers to the 
thermodynamic limit. 

For a fc-regular graph Gt with fc > 3, there is a upper bound for the 
number of spanning trees |23, 123 

NsT{Gk) < ( 4^#^1 {bk^^^^ , (2.2) 



where 



v{Gk)k In k 



|(:(t-2)]i-' 
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With eq. (EH), this then yields [T^] 

ZG, < Mhk) ■ (2.4) 

The construction of the two-dimensional Sierpinski gasket SG2{n) at 
stage n is shown in Fig. ^ At stage n = 0, it is an equilateral triangle; while 
stage n + 1 is obtained by the juxtaposition of three n-stage structures. In 
general, the Sierpinski gaskets SGd can be built in any Euclidean dimension 
d with fractal dimensionality D = \n{d + l)/ln2 For the Sierpinski 

gasket SGd{n), the numbers of edges and vertices are given by 

e(5G,(n))= (^'^+^)(d+l)" = ^(d+ir+i , (2.5) 

v{SGd{n)) = ^^[{d+l)- + l] . (2.6) 

Except the (d+l) outmost vertices which have degree d, all other vertices 
of SGd{n) have degree 2d. In the large n limit, SGd is 2c?-regular. 



A A 

SG2{0) 5G2(1) SG2{2) SG2(3) 

FIG. 1: The first four stages n = 0,1,2,3 of the two-dimensional Sierpinski gasket 
SG2{n). 

The Sierpinski gasket can be generalized, denoted as SGd,b{n)^ by intro- 
ducing the side length h which is an integer larger or equal to two jJO] . The 
generalized Sierpinski gasket at stage n -I- 1 is constructed with h layers of 
stage n hypertetrahedrons. The two-dimensional SG2,b{n) with & = 3 at 
stage n — 1,2 and = 4 at stage n = 1 are illustrated in Fig. |21 The 
ordinary Sierpinski gasket SGd{n) corresponds to the b ~ 2 case, where 
the index b is neglected for simplicity. The Hausdorff dimension for SGd,b 
is given hy D = ln(''+j"^)/ln6 Notice that SGd,b is not /c-regular 

even in the thermodynamic limit. 
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5'G2,3(1) 



SG2,3(2) 



5G2,4(1) 



FIG. 2: The generalized two-dimensional Siorpinski gasket SG'z.bip-) with fe = 3 at stage 
n = 1 , 2 and b = 4 at stage n = 1 . 

III. THE NUMBER OF SPANNING TREES ON SG2(n) 

In this section we derive the number of spanning trees on the two- 
dimensional Sierpinski gasket SG2{n) in detaiL Let us start with the defi- 
nitions of the quantities to be used. 

Definition III.l Consider the two-dimensional Sierpinski gasket SG2{n) 
at stage n. (a) Define f2{n) = NsT{SG2{n)) as the number of spanning 
trees, (b) Define ga2{n), gb2{n), gc2{n) as the number of spanning sub- 
graphs with two trees such that one of the outmost vertices belongs to one 
tree and the other two outmost vertices belong to the other tree as illus- 
trated in Fig. 13 (c) Define h2{n) as the number of spanning subgraphs 
with three trees such that each of the outmost vertices belongs to a different 
tree. 

The quantities /2('t-), 90.2(1^), gb2{n), gc2{n) and h2{n) are illustrated in 
Fig. 13 where only the outmost vertices are shown. It is clear that the 
values ga2 (n) , gb2 (n) , gc2 (n) are the same because of rotation symmetry, 
and we define 52 = go.2{n) = gb2{n) — gc2{n). The initial values at stage 
are /2(0) = 3, (72(0) = 1, /i2(0) — 1. Notice that /i2(n) is the number of 
spanning trees on SG2{n) with the three outmost vertices identified. The 
purpose is to obtain the expression of /2('^) as follows. 

Theorem III.l The number of spanning trees on the two-dimensional 
Sierpinski gasket SG2{n) at stage n is given by 



f2{n) = 2' 



iQ2(")o/32(n)t;72(n) 



(3.1) 



/2(n) gci2(n) ff&2(n) gc2{n) /i2(n) 

FIG. 3: Illustration for the spanning subgraphs /2(n), 902(71), gb2{n), gC2{n) and h2{n). 
The two outmost vertices at the ends of a solid line belong to one tree, while the two 
outmost vertices at the ends of a dot line belong to separated trees. 

where the exponents are 

«2(n) - i(3" - 1) , (3.2) 
/32(n) = ^(3"+i+2n+l) , (3.3) 

72(n) = i(3"-2n-l) . (3.4) 

This theorem can be proved by the following two lemmas. The three 
quantities f2{n)^ Q2{n) and h2{n) satisfy recursion relations which were 
first obtained in ^24J . 

Lemma III. 2 For any non-negative integer n, 

/2(n + l)-6/|(n)52(ri) , (3.5) 

g2{n + 1) = fl{n)h2{n) + 7f2{n)gl{n) , (3.6) 

h2in + 1) = I2f2in)g2in)h2in) + lAgHn) . (3.7) 

Proof The Sierpinski gaskets SG2{n + I) is composed of three SG2{n) 
with three pairs of vertices identified. To obtain the number of spanning 
trees f2{n + 1), one of the SG2{n) should be spanned by two trees. There 
are six possibilities as illustrated in Fig. ^ Therefore, we have 

f2{n + 1) = 2fl{n)[ga2{n) + gb2{n) + gc2{n)] = 6/|(n)g2(n) . (3.8) 

Similarly, 302 (n + 1) for SG2(n + 1) can be obtained with appropriated 
configurations of its three constituting SG2{n) as illustrated in Fig. 
Thus, 



502(71 + 1) = fl{n)h2{n) + 'if2{n)ga\{n) + 2f2{n)ga2{n)gc2{n) 






FIG. 4: Illustration for the expression of /2(n + 1). 



+ 



+ 




+ 



+ 



5:a 



FIG. 5: Illustration for the expression of ga2{n + 1). 



+ 2/2 (71)502 (n)g62(n) 



(3.9) 



With the identity ga2(n) = 562(^1) — gc2{n) — 52 (tt-), eq. (|3.6|) is verified. 

Finally, h2{n + 1) is the summation of appropriated configurations as 
illustrated in Fig. El so that 

h2{n + 1) 

= '^f2{n)h2{n)[ga2{n) + gb2{n) + gc2{n)] 

+2gc2{n)ga2{n)[gc2{n) +502(71.)] + 2ga2{n)gb2{n)[ga2{n) + gb2{n)] 
+2gb2{n)gc2{n)[gb2{n) + gc2{n)] + 2ga2{n)gb2{n)gc2{n) . 

(3.10) 

With the identity 502(71) = 5^2 (t^) = gc2{n) — 52("-), eq. (|3.7|l is verified. 
□ 



x3 + 



x3 + 



x3 



x3 



+ 



+ 



FIG. 6: Illustration for the expression of h2{n + 1). The multiplication for the eight 
configurations on the right-hand-side corresponds to three possible orientations. 
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By the three equations H3.5|l , H3.6|l and H3.7() , /2 (n) , 52 (n) and /12 (f^) can 
be solved. 



Lemma III. 3 For any non-negative integer n, 

/2(n) = 2"^(")3'^^(")5''^(") , (3.11) 

c,2(n) = 2"^(")3'^2(")-"-i5''^(")+" , (3.12) 

/i2(n) = 2"^(")3''^(")~2"~i5''^(")+2" . (3.13) 

where 02 (n), /32(n) and 72(«) are given in eqs. (|3.2() - (|3.4|l 

Proof As 02(0) = 0, /32(0) = 1 and 72(0) = 0, eqs. jTTT|l - PTl|l are 
correct for n ~ 0. By the recursion relation H3.5|l . it is easy to obtain 



where 



/(n+ 1) = 2''"^(")+-'-3''^^2(")^"5'^'''2'^")"'"" 



3a2(n) + 1 = -(3"+i - 1) = a2(n + 1) , 



3/32(n) - n = -(3"+2 + 2n + 3) = ^2(n + 1) 



(3.14) 

(3.15) 
(3.16) 



372(71) +n= -(3"+i - 2n - 3) = 72(71 + 1) 



(3.17) 



The proof for eq. (|3.11|) is completed by induction. Similarly, we have 

g(n + 1) = 2^"=(")+^3^*(")"2""25^''='(")+^"+^ , (3.18) 



h[n + 1) = 23"2(n) + l33/32(n)-3n-3g372(n)+3n-|-2 

which verify eqs. and (fTT^ . □ 



(3.19) 



Although the Laplacian matrix Q{G) for G — S* 6*2(71) can be con- 
structed, it does not look simple to diagonalize. For example, with ap- 
propriate order of the vertices, 



Q(5G2(1)) 



/ 4 -1-1 -1 -1 \ 
-1 4-1-10 -1 

-1 -1 4-1-10 
0-1-12 
-10-10 2 

V -1 -1 2 y 



(3.20) 
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and 

Q(5G2(2)) 
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(3.21) 

Denote the non-zero eigenvalues of Q{SG2{n)) as \{SG2{n))i for 1 < i < 
viSG2{n)) - 1. 



Corollary III.l The product of non-zero eigenvalues of Q{SG2{n)) is 
given by 

v{SG2{n))-l 

n A(5G2(n)), = -(3" + l)/2(n) , (3.22) 

1=1 

where f2{n) is given in Theorem llll. H and the number of vertices v{SG2{n)) 
for SG2{n) is given by eq. (|2.6|l with d — 2. 

By the definition in eq. H2.1() . we have the following corollary. 

Corollary III. 2 The asymptotic growth constant for SG2 is given by 

ZSG2 = ^ + ^ + ^ ^ 

~ 1.048594856... (3.23) 

This is equivalent to eq. (20) of [2^. In passing, we notice that the 
number of spanning trees is the same for a planar graph and its dual, so 
that Theorem llii. II also applies to the dual of SG2{n), denoted as SG^in). 
As SG2 is 4-regular in the large n limit, we have zgQ* = zgQ^ [T^ . 
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IV. THE NUMBER OF SPANNING TREES ON SG2,b{n) WITH 

6 = 3,4 



The method given in the previous section can be appUed to the number 
of spanning trees on SGd,b{n) with larger values of d and b. The num- 
ber of configurations to be considered increases as d and h increase, and 
the recursion relations must be derived individually for each d and h. In 
this section, we consider the generalized two-dimensional Sierpinski gasket 
SG2,b{n) with the number of layers h equal to three and four. For SG2,3{n), 
the numbers of edges and vertices are given by 

e{SG2,3{n)) = 3x6", (4.1) 



v{SG2An)) = ^ , (4.2) 

where the three outmost vertices have degree two. There are (6" — l)/5 
vertices of SG2,3{n) with degree six and 6(6" — l)/5 vertices with degree 
four. Define /2,3(n), g2,3{n), /i2,3('^) as in Definition IIII.ll such that the 
initial values are /2,3(0) = 3, (72,3(0) = 1, ^2.3(0) = 1. By the method 
illustrated in the above section, we obtain following recursion relations for 
any non-negative integer n. 

f2.3{n + 1) = Wl3in)92An)h2An) + 142/1,3(^^)32' 3 (") , (4-3) 



32,3(^ + 1) = 2/4,3 (n) ^3 (n) + 77/, 



2,3Wff2,3W^2,3(»^) 



171/2,3(")52,3(") 



(4.4) 



h2,3in + 1) = Q0flsin)g2Mn)hl^{n) + 564/|,3(n).g3 3(n)/i2,3(ri) 

+ A68f2,3{n)gh{n) . (4.5) 

The figures for these configurations are too many to be shown here. By 
induction as in Lemma IlII. 31 these equations can be solved. 



Theorem IV. 1 The number of spanning trees on the two-dimensional 
Sierpinski gasket 56*2,3 (n) at stage n is given by 

/2 3(n) = 2"^'^'^"^3''^'^'^"^5'*'^'^^"^7''^'^^"^ , (4.6) 
where the exponents are 

a2,3(n) = ^(6" - 1) , (4.7) 
5 
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/32.3W = 4(13x6"-15n+12) , (4.8) 
25 



72,3(n) = ^(3x6"-15n-3) , (4.9) 



'52,3(n) = 4(7x6" + 15"-7) • (4.10) 
25 



52,3 ("■) and h2^j,{n) can also be expressed by these exponents: 

/l2 3(n) — 2"2.3(")3/32,3(«)+2n-l572.3(n)+2ny(52,3(")-2n _ (4.11) 

By the definition in eq. I|2.1|) . we have the following corollary. 

Corollary IV. 1 The asymptotic growth constant for 56*2,3 is given by 
= ^In2+i^ln3+^ln5 + iln7 
~ 1.133231895... (4.12) 

For SG2,i{n), the numbers of edges and vertices are given by 

e{SG2,i{n)) = 3 X 10" , (4.13) 

v[SG2An)) - i^ii^ , (4.14) 

where the three outmost vertices have degree two. There are (10" — l)/3 
vertices of SG2,4{ii) with degree six, and (10" — 1) vertices with degree 
four. Define f2,i{n), 32, 4(h), h2,4{n) as in Definition IIII.ll such that the 
initial values are /2,4(0) = 3, 52,4(0) = 1, /i2,4(0) = 1. We write a computer 
program to obtain following recursion relations for any non-negative integer 
n. 

/2,4(n + l) = 2fl^{n)hliin) + 516flAn)gl^in)hl^{n) 

+ 5856/|4(n)5|4(n)/i2,4(n) + 11354/44(^)32% (n) (4.15) 



92A{n + 1) = 82/2%(n).g2,4(n)/ii4(") + 2786fl^{n)gl^{n)hljn) 
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+ 14480/44(n)g|4N/i2,4H + 13732/2%W32,4W (4-16) 

h2A{n + l) = 20fl^in)hl^in) + 2388fli{n)gl^{n)hl^{n) 

+ 30948/2%(n)5^^(n)/i2_4(n) + 83234/2%(n)5^_4(n)/i2,4(n) 

+ 42210/|,4(n)52%(«) • (4.17) 
By induction as in Lemma [111. 31 these equations can be solved. 

Theorem IV. 2 The number of spanning trees on the two-dimensional 
Sierpinski gasket SG2Aiji) at stage n is given by 

/2,4(n) = 2"^'''(")3''^'^(")5''^'-'("Ul''^'''(")l03"^-*(") , (4.18) 
where the exponents are 

a2A{n) - ^(10" - 1) , (4.19) 

/324N = ^(10" + 2) , (4.20) 

72,4W - ^(10" - 1) , (4.21) 

'52,4(n) = ^(2 X 10" + 9n-2) , (4.22) 

e2,4(«) = ^(10"-9n~l) . (4.23) 

52.4 (f^) and h2^i{n) can also be expressed by these exponents: 
92 iin) — 2"^'''*^"^3''^''''"^~"^5'''^'''''"Ul''^'''''"-'~"103^^'''^"^"''" 

h2 i{n) — 2"2'''*^")3'^=''''*^")~-^5''^''''^"Ul''^''''^"^"^"103^^-'''"''+^" . (4.24) 
By the definition in eq. p.l() . we have the following corollary. 

Corollary IV. 2 The asymptotic growth constant for SG2,4 is given by 
0SG.,4 = iln2+iln3 + ^ln5+iln(41) + ^ln(103) 

~ 1.194401490... (4.25) 
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V. THE NUMBER OF SPANNING TREES ON SGd{n) WITH 

d = 3,4 

In this section, we derive the number of spanning trees on SGd{n) with 
d = 3,4. For the three-dimensional Sierpinski gasket SGsin), we use the 
following definitions. 

Definition V.l Consider the three-dimensional Sierpinski gasket SG3{n) 
at stage n. (a) Define /^{n) = NsT{SG3{n)) as the number of spanning 
trees, (b) Define gsin) as the number of spanning subgraphs with two trees 
such that one of the outmost vertices belongs to one tree and the other three 
outmost vertices belong to the other tree, (c) Define /i3(n) as the number 
of spanning subgraphs with two trees such that two of the outmost vertices 
belong to one tree and the other two outmost vertices belong to the other 
tree, (d) Define (n) as the number of spanning subgraphs with three trees 
such that two of the outmost vertices belong to one tree and the other two 
outmost vertices separately belong to the other trees, (e) Define 93(71) as 
the number of spanning subgraphs with four trees such that each of the 
outmost vertices belongs to a different tree. 

The quantities f^in), gsin), h^{n), pz{n) and 93(71) are illustrated in Fig. 
[7| where only the outmost vertices are shown. There are four equivalent 
(73(71), three equivalent h^{n), and six equivalent p^in). The initial values 
at stage are ^(O) = 16, 33(0) = 3, hsiO) = 1, ^3(0) = 1, 93(0) = 1. 
Notice that 93 (n) is the number of spanning trees on 56*3(71) with the four 
outmost vertices identified. 




/a (71) 53(71) hsin) p:i{n) 173(71) 

FIG. 7: Illustration for the spanning subgraphs fzin), gi(n), h-^in), P3{n) and 53(71). 
The two outmost vertices at the ends of a solid line belong to one tree, while the two 
outmost vertices at the ends of a dot line belong to separated trees. 

In the process of calculation, we find that it is convenient to combine 
(73(71) and h^ln) and define ghzin) = g^(n) + h3{n). We obtain following 
recursion relations for any non-negative integer n. 

hin + 1) = 72/|(7i).g/i3(7i)p3(n) + 56/3(71)5/1^(71) , (5.1) 



5/73(11 + 1) = 6/|(ii)5/i3(n)(73(7i) + 26/|(7i)p2(„) 
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+ l2Qh{n)ghl{n)p:,{n) + 22ght{n) , (5.2) 

P3(n + 1) = &fi{n)pz{n)qz{n) + 14/3(72)3/12(71)93(71) 

+ l2Qh{n)gh:,{n)pl{n)+^&ghl{n)p^{n) , (5.3) 

93(77 + 1) = 144/3(77)5/73(77)^3(77)93(71) + 208/3(r7)p3(n) 

+ 56g/73(7i)q3(n) + 7205/7^(77)^2 („) . (5.4) 

The figures for these configurations are too many to be shown here. With 
the initial values given above, these equations can be solved by induction 
as in Lemma rill.3l 

Theorem V.l The number of spanning trees on the three-dimensional 
Sierpinski gasket SG^{n) at stage 77 is given by 

/3(77) =2"3(")3'^='(") , (5.5) 

where the exponents are 

03(77) =4"+i+ 77 , (5.6) 
/33(77) - i(4"-3n-l) . (5.7) 
5/73(77), ^3(77) and 93(77) can also be expressed by these exponents: 

5/73(77) 2"^(")-"-23/33(n)+« ^ 

^3(77) = 2"-'("^^2"-43/33(n)+2n 

93(77) = 2"^'(")-3«-43/33(n)+3n (5 

By the definition in eq. H2.1fl . we have the following corollary. 



Corollary V.l The asymptotic growth constant for SG^ is given by 

1 
6 

1.569396409... (5.9) 



ZSG3 = 21n2+ -ln3 
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For the four-dimensional Sierpinski gasket SGi{n), we use the following 
definitions. 



Definition V.2 Consider the four-dimensional Sierpinski gasket SGi{n) 
at stage n. (a) Define fi{n) = Nsj'{SGi{n)) as the number of spanning 
trees, (b) Define g4{n) as the number of spanning subgraphs with two 
trees such that two of the outmost vertices belong to one tree and the 
other three outmost vertices belong to the other tree, (c) Define /i4(n) 
as the number of spanning subgraphs with two trees such that one of the 
outmost vertices belong to one tree and the other four outmost vertices 
belong to the other tree, (d) Define P4(n) as the number of spanning 
subgraphs with three trees such that one of the outmost vertices belong 
to one tree, two of the other outmost vertices belong to another tree and 
the rest two outmost vertices belong to the other tree, (e) Define (74 (n) 
as the number of spanning subgraphs with three trees such that three of 
the outmost vertices belong to one tree and the other two outmost vertices 
separately belong to the other trees, (f) Define r4(n) as the number of 
spanning subgraphs with four trees such that two of the outmost vertices 
belong to one tree and the other three outmost vertices separately belong 
to the other trees, (g) Define 54(71) as the number of spanning subgraphs 
with five trees such that each of the outmost vertices belongs to a different 
tree. 

The quantities f 4,(71), gi^n), ft.4(n), Piin), 94(71), r4(n) and 34(71) are 
illustrated in Fig. |H1 where only the outmost vertices are shown. There 
are ten equivalent 34(71), five equivalent /i4(7i), fifteen equivalent ^3(71), ten 
equivalent 94(71) and ten equivalent ^4(71). The initial values at stage 
are /4(0) = 125, 54(0) = 3, h4{Q) = 16, P4(0) = 1, 94(0) = 3, 74(0) = 1, 
54(0) = 1. Notice that 34(77) is the number of spanning trees on 804(71) 
with the five outmost vertices identified. 

We find that it is convenient to reduce the number of variables by defining 
(7/14(71) = 3^4(71) -1-/14(71) aiidpq4(n) = 2^4(77) -1-94(71). We write a computer 
program to obtain following recursion relations for any non-negative integer 
n. 



/4(n+l) 



lU0f^{n)gh4{n)pq4(7i)r4{n) + 520f^(7i)pql(7i) 



+ 1120/4(77)5/71(77)7(77) -h 3600/4(77)5/7^(77)^9^(77) 



+ 1320.9/7^(77)^94(77) , 



(5.10) 



5/74(77 + 1) 



72/4(77)5/14(77)^94(77)54(77) -I- 378/4(77)5/74(77)74(77) 
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ftin) gi{n) h4,{n) P4{n) 




qiin) r4(n) S4(n) 

FIG. 8: Illustration for the spanning subgraphs /4(n), g4{n), h4(n), p4(n), q4{n), r4(n) 
and S4(n). The two outmost vertices at the ends of a solid line belong to one tree, while 
the two outmost vertices at the ends of a dot line belong to separated trees. 

+ 816f^{n)pql{n)rA{n)+56U{n)ghl{n)sA{n) 

+ 3756 f4{n)ghl{n)pq4{n)ri{n) + 2360 f4{n)gh4{n)pql{n) 

+ 688ghlin)r4{n) + 2562ghlin)pql{n) , (5.11) 

pq4{n+l) = 4&fl{n)gh4{n)r4{n)s4{n) +b2fl{n)pql{n)s4{n) 

+ 5Ufi{n)pq4{n)rl{n) + 240f4{n)ghl{n)pq4{n)s4{n) 

+ 1252f4in)ghl{n)rlin) + 4720f4{n)gh4{n)pql{n)r4{n) 

+ 724f4{n)pqtin) + Ughtin)s4{n) 

+ 3A16ghl{n)pq4in)r4{n) + 3104ghl{n)pql{n) , (5.12) 

nin + 1) = 72f^{n)pq4{n)ri{n)s4{n) + 126/|(n)r|(n) 

+ 168/4(n)(?/i4(n)r4(n)s4(n) +360f4{n)gh4{n)pq4{ri)S4{n) 

+ 3756f4{n)gh4{n)pq4{n)rl{n) + 2360 f4{n)pql{n)r4{n) 

+ 26Aghl{n)pq4in)s4in) + 1376ghl{n)rl{n) 

+ 7686ghl{n)pql{n)r4{n) + 2328gh4{n)pql{n) , (5.13) 

S4(n + 1) = 2880f4{n)gh4{n)pq4{n)r4{n)s4{n) + 504:0f4{n)gh4{n}rl{n} 
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+ lOAO Mn)pql{n)s4n) + 16520 U{n)pql{n)rl{n) 

) + 3600ghl{n)pql{ 
+ 37560ghl{n)pq4,{n)rl{n) + A7200ghi{n)pql{n)r4,{n) 

+ 43Mpql{n) . (5.14) 

With the initial values given above, these equations can be solved by 
induction as in Lemma IlII. 31 

Theorem V.2 The number of spanning trees on the four-dimensional Sier- 
pinski gasket SGi{n) at stage n is given by 

fi{n) = 2"''(")5''*(")7'^''(") , (5.15) 
where the exponents are 

«4(r^) = ^(5" ~ 1) , (5.16) 

/34(n) = ^(5"+i+47^ + 3) , (5.17) 

o 

74(n) = ^(5"-4n-l) . (5.18) 

(7/14(71), pqi{n), r^ln) and 54(71) can also be expressed by these exponents: 
gh^n) = 2"''(")5'^*("^-"-^7'^''(")+" , 

pqi{n) = 2"*(")5'^*("'^^"^^7'''^("^+^" , 

r4(7i) = 2"^*^"'5'^*'^"'^'^"^'^7''*("^+^"' 

54(71) = 2"''(")5''''(")-'*"-37''''(")+''" . (5.19) 
By the definition in eq. H2.1|) . we have the following corollary. 

Corollary V.2 The asymptotic growth constant for SG4 is given by 
ZSG, = ^In2 + ^ln5 + |jln7 

~ 1.914853265... (5.20) 
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VI. THE NUMBER OF SPANNING TREES ON SGd{n) FOR 

GENERAL d 

With the above theorems that the numbers of spanning trees on SG2 (n) , 
S* 6*3(71) and SGi(n) have only simple factors, we observe the following 
conjecture for the number of spanning trees fdin) = NsriSGdin)) on 
SGd{n) for general d. 



Conjecture VI. 1 The number of spanning trees on the d-dimcnsional 
Sierpinski gasket SGd{n) at stage n is given by 

fd{n) = 2"''(")(d + l)'3''(")(d + 3)'^'^(") , (6.1) 

where the exponents are 

a4n) = ^[id+ir-l] , (6.2) 



Mn) - ^ [{d + 1)"+^ + dn + d-l], (6.3) 



^'ii^)-^[id+^r-dn-l] . (6.4) 



Notice that these exponents are positive integers when d is a positive 
integer and n is a non-negative integer. It is interesting that the recursion 
relations l|n^ - l|lT7jl . ifOjl - lfCTjl and lfOn|l - ((OI|l are very 

complicated but the solutions are very simple. We do not know how to 
explain it. There may be a better method to solve this problem. 

SGd{0) at state n = is a complete graph with {d + 1) vertices, each of 
which is adjacent to all of the other vertices. Setting n — into Conjecture 
IVI.ll is consistent with the expectation that the number of spanning trees 
on the d-dimensional Sierpinski gasket SGd{0) at stage zero is given by 

fd{o) = {d+ir-\ 

With the number of vertices for SGd{n) given in eq. H2.6|l . Conjecture 
IVI. II leads to following result. 



Conjecture VI. 2 The asymptotic growth constant for SGd{n) is given by 

d- 
d(d 



zsG, = rfM \ [dln2 + id+l) ln{d + 1) + ln(rf + 3)] . (6.5) 



18 



TABLE I: Numerical values of zsoa s-nd ^SGa- The last digits given are rounded 
off. 



d 


D 


k 






2 


1.585 


4 


1.048594857 


0.862051042 


3 


2 


6 


1.569396409 


0.928042816 


4 


2.322 


8 


1.914853265 


0.953722370 


5 


2.585 


10 


2.172764568 


0.966999152 


6 


2.807 


12 


2.378271274 


0.974945363 


7 


3 


14 


2.548944395 


0.980152960 


8 


3.170 


16 


2.694814686 


0.983785240 


9 


3.322 


18 


2.822140640 


0.986437211 


10 


3.459 


20 


2.935085659 


0.988442577 



Co nj ect u re IVLTI and IVO reduce to Theorems ITTlTI lyU and Corol- 
laries |^OJ ^3 ^3 respectively when d is set to two, three and four. 

By eq. (|2.4I) . zsOd has the upper bound ln52rf. It is of interest to see 
how close the value zsCd to this bound. For this purpose, we define the 
ratio 

rsG, = ff^ , (6.6) 

and list the first few values of zsOd JsGj h^ Table Our results 
agree with the observation made in Ref. [ill ^^^^ ^SGd increases as the 
degree k = 2d increases. Compared with the values zr, for d-dimensional 
hypercubic lattice £d which also has fc = 2c? in 0, [sy , zsg^ < for 
all d > 2 indicates that SGd is less densely connected than Cd- The ratio 
rsGd increases with dimension and approaches one, but the convergence is 
not as fast as those for the regular lattices given in 0, . 
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